Magnetic quantum oscillations in doped antiferromagnetic insulators 
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Energy spectrum of electrons (holes) doped into a two-dimensional antiferromagnetic insulator is 
quantized in an external magnetic field of arbitrary direction. A peculiar dependence of de Haas- van 
Alphen (dHvA) or Shubnikov-de Haas (SdH) magneto-oscillation amplitudes on the azimuthal in- 
plane angle from the magnetization direction and on the polar angle from the out-of-plane direction 
is found, which can be used as a sensitive probe of the antiferromagnetic order in doped Mott- 
Hubbard, spin-density wave (SDW), and conventional band-structure insulators. 
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Quantum oscillations of magnetization and resistivity 
with the magnetic field are of a great experimental and 
theoretical value providing reliable and detailed Fermi- 
surfaces [H, H, [1| . Specifically interest in dHvA and SdH 
effects in almost two-dimensional (2D) Fermi-liquids has 
recently gone through a vigorous revival due to experi- 
mental discoveries of magneto-oscillations in a few high- 
temperature cuprate superconductors [1, d, d, 0] • Their 
description in the framework of the standard theory for 
a metal [l[ has led to a small electron-like Fermi-surface 
area of a few percent of the first Brillouin zone and 
to a surprisingly low Fermi energy of about the room 
temperature [6), |7|] , somewhat inconsistent with the first- 
principle (LDA) band structures and angle-resolved pho- 
toemission (ARPES) spectra of cuprates Q. The oscilla- 
tions have been observed in the superconducting (vortex) 
state well below the upper critical field raising a doubt 
concerning their normal state origin While a better 
understanding of dHvA/SdH effects in doped antiferro- 
magnetic insulators is generally important, it becomes 
particularly vital for building an adequate theory of high- 
temperature superconductivity since parent cuprates are 
antiferromagnets. 

Here, using a tight-binding Hamiltonian we quantize 
the energy spectrum of electrons or holes moving on 
the anti-ferromagnetic (AF) background in a two dimen- 
sional lattice. We find a peculiar dependence of the 
magneto-oscillation amplitudes on the magnetic field di- 
rection, which could serve as a sensitive probe of the 
antiferromagnetic order in doped insulators. 

The mean-field tight-binding Hamiltonian of carriers 
doped into the bipartite antiferromagnetic square lattice 
in the external magnetic field, B, is written as 

H = ^^Sii^^^ala^Gi + fiBBalerai) + tii>alai 
a' 

+ YUjb]a, + H.c, (1) 

where a\ — (a||, a\t^) and b\ — [b\.^, b\^) create the carrier 
on sites "i" and "j" of sublattices A and B, respectively, 
with the spin s ^],[, A is the carrier spin-lattice spin 



FIG. 1: The azimuthal in-plane angle, $, from the magneti- 
zation direction, n, and the polar angle, 0, of the magnetic 
field B from the out-of-plane direction. 



exchange energy (the antiferromagnetic gap), Sa' is the 
Kroneker symbol, tu' , tjji and tij are the hopping inte- 
grals, and cr = {cr^;, Uy, Cz} are the Pauli matrices. 

Fourier transforming the operators from Wannicr (site) 
to Bloch (momentum, k) representation and assum- 
ing translational invariance the carrier energy spectrum, 
-B(k), is found by diagonalizing 4x4 matrix: 



4 - Actq - ^b5[|0-2 liBB±ao + ti^ax 



(2) 



where ij^ = I],,' *3xp(i[k • (j' - j)]) is the hopping 



energy within one sublattice, tk = ^j^- exp(z[k • (i — 
j)]) is the inter-sublattice hopping energy. This matrix 
corresponds to the choice of the 4-dimensional vector in 
the spin and sublattice space at fixed k. Here B± and 
-B|| are transverse and longitudinal components of the 
magnetic field with respect to the lattice magnetization 
n, Fig.l, (Tq is the identity matrix, and hb is the Bohr 
magneton. There are two electron and two hole bands 
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dispersed as 



Bi 



1/2 



(3) 



in the lowest order with respect to the field, fisB <C A. 
They are split by the external magnetic field into two 
subbands each with anysotropic g-factor, g — 2 [cos^ (</>) + 
sin^((/))t^/(A^ + tU] depending on the angle (j) between 
the field and the magnetization, Fig.l. 

The anisotropic g-factor differs significantly from the 
free-electron = 2 near the extremum points of the va- 
lence/conductance bands, where <C A^. According 
to general principles of quantum mechanics deviations of 
the ^-factor from its classical value are related to spin- 
orbit interaction. The spin-orbit interaction is not in- 
cluded explicitly to the Hamiltonian Eqs.(l,2). Basically 
the difference originates from the spin-orbit interaction 
pinning the lattice magnetization along a crystal lattice 
direction and present in the Hamiltonian implicitly. At 
a relatively low doping with the Fermi energy, Ep near 
the top (bottom) of the valence (conduction) band, one 
can expand Eq.Q in powers of t/A, 



E{k) 



2mr 



2m,, 



± Ms 



Bf,+j^(k)B 
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(4) 



Here m,-^ = Aa'^{2t^/A - t)/h^, rUy^ = Aa^t /h^ are 
components of the effective mass tensor, a is the lattice 
constant, t and t > are nearest and nearest next neigh- 
bor hopping integrals, respectively, and the coefhcient 
7(k) is small as 7(k) = 2V^tK/l\ - {Ep/AY''^ < 1. 

are deviations of the wave vector perpendic- 
ular and parallel to the antiferromagnetic Brillouin zone 
boundary, respectively, and the energy of the extremum 
point is taken as zero. 

The anisotropic (j-factor in doped antiferromagnetic in- 
sulators was originally derived in a weak-coupling nesting 
model fio']. Actually the effective mass approximation, 
Eq. Q , can be also derived phenomenologically using the 
symmetry arguments (llj . The non- unitary group of the 
antiferromagnetic lattice is G = {£)4,i?T}, here de- 
scribes all rotations which remain the system invariant. 
Translations T by a lattice period transform from one 
sublattice to another changing the sign of the magneti- 
zation. Hence, these translations are multiplied by the 
time inversion operator R. Following Brazovskii and 
Lukyanchuk JTj one can construct the Hamiltonian of 
the required symmetry as 
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H 
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I a'(Toa~ 



V 



[(B • n)a^(n • cT)a + 7(k)(B x n)a){n x cr)a](5) 

with the electron (hole) energy spectrum Eq.Q. Here 
n is the magnetization unit-vector, and a) = {a\^a}^),a 
are creation and annihilation operators, respectively, for 
the spinor describing the hole (electron) band. The co- 
efficient 7(k) is an odd function of kx, which is zero at 



the antiferromagnetic Brillouin zone boundary, so that 
7(k) = ^kxi where 7 does not depend on k. The cou- 
pling to the magnetic field in this Hamiltonian is ob- 
tained noticing that the transformation k ^ k + Q with 
Q = ira/a? is equivalent to the rotation in the spinor 
space described by the matrix n cr [ll[ (here a = {a, a}). 
Direct comparison of the eigenvalues of the Hamiltonian 
Eq.(5) with the spectrum Eq.(3) yields 7 = 2^/21/ A. Im- 
portantly, the symmetry arguments are applied beyond 
the mean- field approximation, Eq.(IT]), so that spin fluc- 
tuations just renormalize the effective mass tensor and 
other coefficients in Eq.(|3]). 

The orbital quantization of the spectrum, Eq.Q, is 
readily obtained via the Peierls substitution [IJl, k =J> 
— iV + eA with the vector potential A(r) in Eq.(2). In 
the lowest order with respect to EpjA. we can use the 
effective mass approximation, E g .(141). which yields the 
conventional Fock-Landau levels [ij, [lj| split by the lon- 
gitudinal field as 



£;„ = cos(e)|(n 1/2) ± 



(6) 



where w — eB j (raxmy)^ 1"^ is the cyclotron frequency, n = 
0, 1,2..., and 8 is the polar angle between the magnetic 
field and the out-of-plane direction, Fig.l. 

Now the oscillating part of the magnetization, M, is 
calculated following the standard route by applying the 
Poisson summation [l|: 



~ . 2TTrF 

^ B 



mxTny)^/"^ tan(6) cos(<I') 



Here 

Mr = Ar(e)cOS 

is the amplitude of r-harmonic with 
^ i^huj cos(0) ) ^ 



(7) 



(8) 



27rrr 
fku cos(0) 
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F — {mxrayY/'^Ep/ehcos{Q) is the fundamental fre- 
quency of oscillations, Rt{z) = z/ sinh(2:) and Rd{z) = 
exp(— z) are conventional temperature and Dingle reduc- 
tion factors, r is the scattering rate, m^ is the free elec- 
tron mass, d is the inter-plane distance, and 4> is the az- 
imuthal in-plane angle from the magnetization direction, 
Fig.l. Both angles 8 and $ in Eq.® are changing in the 
interval 0^8,$^ 7r/2. Three-dimensional corrections 
to the energy spectrum_can be accounted for by the ad- 
ditional Yamaji factor [isj, Ry = Jo [47rrtj_//ia;cos(8)] 
in Eq.Q, where Jo{x) is the zero-order Bessel function, 
t± — t±jQ{kFdtan{Q)), t± is the out-of-planc hopping 
integral, and hkp is the Fermi momentum. 

As follows from Eq.® the essential anisotropy of the 
(/-factor causes a strong dependence of the oscillation am- 
plitude on the azimuthal in-plane angle of the field from 
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FIG. 2: dHvA first-harmonic amplitude as a function of the 
azimuthal in-plane angle from the magnetization direction, <1> 
and the polar angle from the out-of-plane direction, O, in a 
layered antiferromagnet (upper panel) compared with the first 
harmonic amplitude in a nonmagnetic layered metal (lower 
panel) at T = F = 0, and {mxmyY^'^ = me- 



the magnetization direction. Fig. 2a, which is absent in 
ordinary non-magnetic layered metals. Fig. 2b, where the 
magnetization amplitudes are found as 



m;'^'' = A^(e)cos 



7rr(TOj.TOj^)"'^/^ 
me cos(9) 



(10) 



The novel dependence on $ and 8, Eq.(l8]), is extremely 
pronounced at low temperatures (compare Fig. 2 (upper 
panel) and Fig. 2 (lower panel)), as also shown in Fig. 3 
for some fixed azimuthal angles. 

One can readily generalize our results to any shape 
of the Fermi surface, and calculate corrections to am- 
plitudes and fundamental frequencies of higher order in 
Ep/A and in the magnetic field by applying the Lifshits- 
Kosevich quasi-classical approximation [l6i |. Within the 
approximation dHvA frequencies F± are determined by 



FIG. 3: dHvA first-harmonic amplitude as a function of 
the polar angle, O, for two different azimuthal in-plane an- 
gles, $, in a layered antiferromagnet compared with the 
independent first harmonic amplitude in a nonmagnetic lay- 
ered metal at T = 0, 2ttV = hcj, 'ykp = 0.1, and {rnxmy)^^'^ = 
nie. 








FIG. 4: dHvA first-harmonic amplitude as a function of the 
polar angle, O, in a disordered antiferromagnet (solid line) 
compared with the first harmonic amplitude in a nonmagnetic 
layered metal (dotted line) at T = F = 0, and {rnxmyY^^ = 



the extremal cross-section areas, S'lf* of two spin-split 
electron (or hole) Fermi surfaces, F± = hS^'^/2TTe. Fol- 
lowing Ref . [l^ one can expand the extremal cross-section 
area in powers of the magnetic field, so that F± = 
F ± aB + (3B^ ± eB^ . Here the second term describes 
the Zceman splitting of the bands with the anysotropic 
g- factor. It does not shift the frequency but affects the 
amplitude. The third term describes a small shift of the 
fundamental frequency, F, depending on the magnetic 
field. The last term describes a small field-dependent 
correction to the (^-factor. For example, when the field is 
perpendicular to the magnetization, Bh = 0, and the 
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effective mass approximation is applied near X-point, 
(7r/2a, 7r/2a), of the antiferromagnetic Brillouin zone, 
one finds F = hkjp/2e, a — {nixmyy^'^jkp /nrrie, /3 = 

mx^'^ml^'^j'^e/Smlh, and e = mt^'^my^'^j^e'^ /12'Kmlh'^kF 
with fikp = [2{mxmyy^'^ EpY^^ . For an arbitrary field 
direction one obtains, using Eq.(j4]) with 7(k) — ^kx <C 1, 



2r [rrixmyicos^icj)) + -^'^k'], sin^((^))] ^'"^ E[K[(j))\ 



cos 



; cos(9) 



(11) 



Here 4> is the angle between the magnetic field and the 
magnetization, Fig.l, E{k) is the elliptic integral of the 
second kind, and 



sin^(0) 



cos^ {(j)) + 7^/c|n sin (0) 



1/2 



Taking 7 = in Eg. pTjl one obtains Eq.® since 
cos2(^) = sin2(e)cos2($) and E[Q] = 7r/2, E[l\ = 
1. The finite transverse spin-susceptibility, oc jkp = 
(2i?i?/A)^/^, only slightly blurs the strong ^-dependence 
of the amplitudes. Fig. 3, if Ep/A <C 1. For example, 
when the field is rotated in the plane perpendicular to 
the magnetization axis n we have 



= Ar{Q) cos 



2r{2mxmyEp / IS) 



1/2 



; cos(9) 



(12) 



with a small transverse ^-factor. On the other hand if the 
magnetic field is rotating in the (z, n) plane, the angular 
dependence is quite different. 



Mr = A(0)C0S 



■Kr{mx'myY^''^ tan(O) 



(13) 



as in Eq.® with $ = 0. 

Real antiferromagnetic solids, like cuprates, could be 
disordered or twinned, so that the magnetization direc- 
tion n within the plane is random. Nevertheless the de- 
pendence of dHvA amplitudes on the polar angle, 0, re- 
mains rather unconventional. Indeed averaging Eq.® 
over all directions of $ from zero to 7r/2 yields 



{Mr)^Ar{Q)Jo 



T:r{raxray)^/'^ tan(0) 



(14) 



which is distinguishably different from the amplitudes in 
a nonmagnetic metal, Eq. (fTO|) . Fig. 4. There are known 
relations between oscillations in transport and thermody- 
namic quantities [l|, at least in nonmagnetic substances. 
Relying on them, we expect the similar nontrivial angle 
dependences also in the SdH magnetooscillations. 

In summary, we have derived the energy spectrum of 
electrons (holes) doped into a two-dimensional antiferro- 
magnetic insulator in terms of all-neighbours hopping in- 
tegrals of nonmagnetic lattice, Eq.®, and quantized it in 
the external magnetic field of arbitrary direction. The pe- 
culiar dependence of dHvA/SdH magneto-oscillation am- 
plitudes on the azimuthal in-plane angle from the magne- 
tization direction and on the polar angle from the out-of- 
plane direction is found, which could be instrumental as 
a sensitive probe of the antiferromagnetic order in doped 
Mott-Hubbard, spin-density wave (SDW), and conven- 
tional band-structure insulators. 
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